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We propose a very simple experimental setup to measure, 
via photon counting, the overlap of the Wigner functions char- 
acterizing two single mode light beams. We show that this 
scheme can be applied to determine directly the phase space 
quasiprobability distribution of the single mode field and in 
a certain limit the Wigner function can be measured without 
use of tomographic reconstruction algorithms. The deleteri- 
ous effects of non-unit photodetector efficiency are analyzed. 

PACS number(s): 42.50.Dv, 03.65. Bz 



The concept of phase space, fundamental in classical 
mechanics, remains useful when passing to quantum the- 
ory. It is well known that quantum states can be fully 
characterized by the Wigner function denned in the phase 
space and that the quantum expectation values can be 
represented as statistical averages of the phase space vari- 
ables [Q. In quantum optics this formalism provides a 
convenient framework for discussion of many topics 0] . 

An exciting problem explored recently is the experi- 
mental determination of the Wigner function of a single 
light mode. It was first shown theoretically by Vogel and 
Risken ]3| that quadrature amplitude distributions mea- 
sured in homodyne detection provide enough data to per- 
form a complete reconstruction of the Wigner function. 
This method, called quantum tomography, was success- 
fully realized in a series of experiments [|J . The quantity 
recorded in the quantum tomography experiments was 
the statistics of the count difference of photodetectors 
facing the signal field superimposed on a strong coherent 
field. The Wigner function was reconstructed from these 
data using numerical algorithms of the inverse Radon 
transform. There are also other known methods for 
the complete experimental characterization of the single 
mode state: heterodyne || and double homodyne 
detection, where the so-called Q function of the signal 
field, which is a smoothed Wigner function, is measured. 

In this Letter we propose an alternative scheme for 
the phase space measurement. Given two single mode 
light beams we present an extremely simple experimen- 
tal setup to measure the overlap of the Wigner functions 
characterizing these field. Moreover we show that the 
Wigner functions can be relatively rescaled by an arbi- 
trary positive factor. The measurement is performed by 
optical means and only a trivial arithmetic operation has 
to be done on the data recorded from the photodetector. 

Before we discuss the proposed measurement, we re- 
view briefly, using a Heisenberg-picture-type of approach, 



the properties of the phase space distributions used in 
this Letter. The Wigner function or the Q function 
are examples of more general s-parameterized quasiprob- 
ability distributions H^(a; s) in the complex a phase 
space [||, with s=0 corresponding to the Wigner function 
W^a) = W(a;0). These distributions can be written as 
expectation values of the following normally ordered op- 
erator: 



U(a;s) 
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where a) and a are the single mode photon creation and 
annihilation operators. In the proposed measurement 
scheme we utilize the simple fact that the single mode 
Wigner function at the origin of the complex phase space 
W(0) can be directly computed from the distribution of 
counts measured by a photodetector facing this mode. 
Indeed for a = and s — the operator ([[]) can be 
written in the following equivalent forms: 
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where we have used the normally ordered operator rep- 
resentation of the n photon number projection operator. 
The quantum expectation value of this operator gives: 



W(0) = (U(0; 0)> 
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where the value p n appearing in this expansion is just the 
probability of counting n photons by an ideal photodetec- 
tor. Thus the photoncount statistics allows to calculate 
the Wigner function at the origin of the phase space. It 
has been pointed out by Royer |h]] that the ability to 
measure W(0) allows us to scan the complete Wigner 
function by shifting the system or equivalently the frame 
of reference in the phase space. Although in principle 
this shifting of the reference frame might be realized in 
quantum optics, its experimental realization would en- 
counter difficulties. Therefore we will utilize the results 
of the measurement of W(0) in an alternative way, pro- 
viding a much more feasible experimental scheme. The 
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proposed scheme is very general and we believe that it 
may find many various applications and generalizations. 
In the present Letter we will discuss in detail one specific 
example motivated by the quantum tomography experi- 
ments. 

The proposed setup is presented in Fig. |l|. We will 
take the detected field to be a superposition of two single 
mode fields, which we will call the signal and the probe, 
and denote their annihilation operators by as and dp 
respectively. Such a combination can be easily realized 
by means of a beamsplitter. As it is well known 
the action of the beamsplitter is described by an SU(2) 
transformation between the annihilation operators of the 
incoming and outgoing modes. Since the phase shifts 
appearing in this transformation can be eliminated by 
an appropriate redefinition of the modes, the annihilation 
operator of the outgoing mode falling onto the detector 
face a ou t can be assumed to be a combination 

a ou t = VTa s - VT- Tap (4) 

where T is the beamsplitter power transmissivity The 
Wigner function of the outgoing mode at the phase space 
origin is given in terms of the incoming modes by the 
expectation value of 

#out(0;0) = - : exp (^a^Oout) : 

= - : exp f-2T(at - — T)/Ta p ) 

x {as - y/(l - T)/Ta P j) : . (5) 

This simple relation provides an interesting link between 
the detected quantity and the S mode. Let us for ex- 
ample consider a simple case when the probe field is a 
coherent state ap\a) = a\a) uncorrelated with the sig- 
nal mode. Performing the quantum average over the P 
mode in Eq. ([|) is straightforward due to the normal or- 
dering of the operators. Taking the expectation value 
over the signal and using the definition (|l|) we obtain 
that the Wigner function (0) for the outgoing mode is 
proportional to an s = 1 — I /T ordered quasidistribution 
function of the S mode: 

Wout(o) = ^w s (J^-<*i - l ~^r) ■ (6) 

Thus our setup delivers directly the value of the signal 
quasidistribution function at the phase space point de- 
pendent on the amplitude and the phase of the probe 
coherent state. Since both these parameters can be con- 
trolled experimentally without difficulties, we may simply 
scan the phase space by changing the amplitude and the 
phase of the probe field and thus determine the complete 
quasidistribution function. Eq. (^) shows that its order- 
ing depends on the beamsplitter transmissivity and for T 
near one approaches zero, which means that the scanned 



quasidistribution is close to the Wigner function of the 
signal field. 

We shall now generalize Eq. (||) for arbitrary state of 
the probe mode by considering the disentanglement of 
the S and P modes. In order to achieve this we use 
the following Gaussian integration of normally ordered 
operators of the S and P modes: 

£>out(0;0) 

= ^ / & 2 f3 : exp (-2(Vr/T - a P ){VT0 - a P )) : 

x : exp (-2(VT^Tf3* - a^)(\/l — T/3 - as)) : ■ 

(7) 

Under the assumption that the S and P modes are un- 
correlated this disentanglement yields the following ex- 
pression for the quantity detected by our setup: 

W O ut(0) = J d 2 f3Ws{VT^T{3)Wp(VTf3) 

= J d 2 pW s (P)W P (VT/(l-T)p). (8) 

This formula establishes the connection between the pho- 
ton statistics of the outgoing mode and the Wigner func- 
tions of the S and P modes. It reflects the fundamental 
advantages of our setup in the direct probing of the phase 
space of the light field. In the case when the beamsplit- 
ter splits the light equally, we have y/T/(l — T) = 1 
and Wout(O) is simply a doubled overlap of the signal 
and probe Wigner functions. In the general case the 
phase space parameterization of the probe Wigner func- 
tion is rescaled by the factor y/Tj (1 — T) which can take 
an arbitrary positive value depending on the beamsplit- 
ter transmissivity. This rescaling causes an effective de- 
crease of the probe width by the factor y/ {1 — T)/T in all 
the quadratures simultaneously. This contrasts with the 
unbalanced double homodyne detection scheme 
where the resolution of the phase space probing along 
one quadrature can be improved only at the cost of de- 
teriorating the resolution along the perpendicular direc- 
tion. For the coherent probe field |a) the P mode Wigner 
function is of the form 

W P (fi) = - exp (-2|/3 -a\ 2 ). (9) 

An easy computation shows that in this case Eq. (||) 
indeed reduces to Eq. (^). The rescaling of the probe 
Wigner function in the convolution (ra) has another con- 
sequence. When T tends to one, the factor multiplying 
the probe amplitude a becomes very small and to scan 
the interesting region of the signal phase space one has 
to use a probe field of large intensity. 

One of the advantages of balanced homodyne detection 
used in quantum tomography experiments is the cancel- 
lation of the excess noise of the reference field. In our 
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scheme this noise deteriorates the resolution with which 
the signal phase space can be probed. This can become 
important in the limit T — > 1, where strong probe fields 
have to be used. The influence of the excess noise can 
be simply estimated assuming a Gaussian thermal noise 
described by the following Wigner function 



?r(2n+ 1) 
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(10) 



where n is the mean number of thermal photons in the 
beam. An easy calculation shows that for such a noisy 
probe field 



Wout(0) = -W s 
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In particular, when n>l grows linearly with the probe 
intensity |a| 2 , the excess noise imposes a restriction on 
the highest ordering of the signal quasidistribution func- 
tion measured at a given point. 

The proposed setup is an optical realization of a model 
scheme of quantum measurement |^2| , |l3| , where in addi- 
tion to the system a filter device — a "quantum ruler" — 
is introduced and the measured phase space probability 
distribution is the convolution of the system and filter 
Wigner functions. Our scheme is more general, since the 
Wigner function of the filter can be rescaled by an ar- 
bitrary factor. Consequently the rescaled probe Wigner 
function does not have to obey the Heisenberg uncer- 
tainty principle and may even approach the shape of a 
delta function, which leads to the direct measurement of 
the Wigner function. In contrast to quantum tomogra- 
phy no sophisticated computer processing of the experi- 
mental data is necessary. The quantity measured in the 
experiment is proportional to the quasiprobability distri- 
bution at the phase space point depending only on the 
amplitude and phase of the probe state. 

In the remaining discussion we will introduce two gen- 
eralizations. First we will make our considerations more 
realistic by taking into account the imperfectness of the 
photodetector. When the detector efficiency is 77, the 
probability of counting n photons is given by the expec- 
tation value of : exp(-r)al ut a out )(r]al ut a ou t) n / nl :. The 
second extension is the substitution of the factor (—1)™ 
in Eq. (||) by — (s + 1)™/ (s — 1)" +1 , where s is a real pa- 
rameter. The origin and role of the parameters 77 and s is 
different: 77 describes experimental limitations, while s is 
an artificial number introduced in the numerical process- 
ing of the measured data. With these two parameters we 
obtain the following simple generalization of the formula 
(H), when expressed in terms of the S and P modes 
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The third line of this equation suggests that the parame- 
ter s can be used to compensate the imperfectness of the 
photodetector. Indeed if we selected s = 1 — 77, we would 
determine the expectation value of : exp(— 2al ut a ou t) : 
regardless of the detector efficiency. But in this case the 
factor multiplying the probability of counting n photons 
is (l — 2/r7) n anc l it s magnitude diverges to infinity with 
n — > 00. This is not important from a theoretical point 
of view, since we have shown that the series converges to 
the expectation value of a well behaved operator. How- 
ever it becomes crucial in the processing of the measured 
probability distribution, which is influenced by experi- 
mental errors and statistical fluctuations, and hence need 
not tend to zero sufficiently quickly to assure the conver- 
gence of the complete series. This reasoning might be op- 
posed since the experimental sample of the photodetector 
counts is finite and thus the counts distribution is zero 
above certain photon number. Nevertheless the problem 
still exists, since the increasing factor in the generaliza- 
tion of the sum (^|) causes that an important contribution 
comes from the "tail" of the experimental counts distri- 
bution, which has usually has a very poor statistics, and 
consequently the final result has a huge statistical error 
The simplest way to avoid all these problems is to 
assume the factors multiplying the counts statistics to be 
bounded, which is equivalent to the condition s < 0. 

An easy computation shows that if the coherent state 
(0) is employed as a probe, the expectation value of the 
generalized operator [/^(O, s) is again given by the qua- 
sidistribution function of the signal mode: 
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Let us now analyze the ordering of this function. It has 
been pointed out by Leonhardt and Paul pBJ that al- 
though from a theoretical point of view an arbitrarily 
ordered distribution contains the complete characteriza- 
tion of the quantum state, experimental errors make it 
difficult to compute higher ordered distributions from the 
measured one. Thus what is interesting is the highest or- 
dering achievable in our scheme. Analysis of the role 
of the parameter s is the simplest, since the greater its 
value, the higher is the ordering obtained. Because it is 
restricted in its range to nonpositive values, it should be 
consequently set to zero. Thus we are left with two pa- 
rameters: 77 and T. It is easy to check that for fixed 77 
the highest ordering is still achieved when T — > 1 , but its 
limit value is now —(l—r))/r). Under the assumption that 
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r\ and T are close to one, the ordering of the measured 
distribution is effectively equal to this limiting value if 
the difference 1 — T is much smaller than 1 — 77. For cur- 
rently used photodetectors, this condition can be realized 
experimentally. Thus the highest ordering achievable in 
our scheme is effectively determined by the photodctcc- 
tor efficiency and is equal to — (f —T])/r). ft is noteworthy 
that this is exactly equal to the ordering of the distribu- 
tion reconstructed tomographically from data measured 
in the homodyne detection with imperfect detectors |L6| . 

The measurement of the quasiprobability distribution 
does not exhaust possible applications of the proposed 
setup. Since the probe field may be in arbitrary state, the 
variety of information on the quantum state which can 
be retrieved using this scheme is potentially very large. 
Another interesting extension of the presented work is its 
generalization to the multimode case. 

We wish to acknowledge useful discussions with J. 
Mostowski and Cz. Radzewicz. We thank P. L. Knight 
for numerous comments about the final version of the 
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the Polish KBN grant. 

Note added. After this Letter was submitted, Ref. 
was brought to the authors' attention, ft discusses 
an analogous measurement scheme with a coherent state 
used as a probe. 
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FIG. 1. Experimental setup discussed in the Letter. BS 
denotes the beamsplitter, PD is the photodetector, and the 
annihilation operators of the modes are indicated. 
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